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Abstract 

I report on some theoretical progress concerning the calculation of gluon bremsstrahlung for 
very high energy particles crossing a weakly-coupled quark-gluon plasma, (i) I advertise that 
two of the several formalisms used to study this problem, the BDMPS-Zakharov formalism and 
the AMY formalism (the latter used only for infinite, uniform media), can be made equivalent 
when appropriately formulated, (ii) A standard technique to simplify calculations is to expand 
in inverse powers of logarithms ln(£'/r). I give an example where such expansions are found to 
work well for co/T > 10 where a> is the bremsstrahlung gluon energy, (iii) Finally, I report on 
perturbative calculations of q. 



1. Equivalence of BDMPS-Zakharov and AMY 

The calculation of gluon bremsstrahlung from high-energy particles traversing a quark-gluon 
plasma is complicated by the Landau-PomeranchukMigdal (LPM) effect: at high energy, the for- 
mation length for bremsstrahlung (the distance over which gluon or photon emission is coherent) 
becomes longer than the mean free path between scatterings. As a result, bremsstrahlung from 
successive scatterings cannot be treated as independent. For QED, Migdal solved this problem 
in 1956. The basic formalism for treating the LPM effect in QCD was worked out by Baier, 
Dokshitzer, Mueller, Peigne and Schiff (BDMPS) and Zakharov in 1996-1998 fl]. There have 
been a number of variations on attacking the problem since then, but the one I want to focus on 
here is that of myself, Moore, and Yaffe (AMY) [2]. The AMY formalism has been used for the 
calculation of transport coefficients, such as shear and bulk viscosity, to leading-order in a s [3, 4]. 
The AMY formalism was derived independently because, at the time, we did not understand the 
BDMPS-Zakharov formalism and could not see how to overcome some of its limitations. 

One of the limitations of the AMY formalism was that it was formulated only for the case 
of a uniform, infinite-volume system, which means systems which are approximately uniform 
over one formation length (and do not change over the corresponding formation time). BDMPS- 
Zakharov, in contrast, could handle the non-uniform, time-dependent case as well, such as a 
system whose size is smaller than the formation length. On the downside, BDMPS-Zakharov 
implicitly treated plasma particles as static scatterers. Furthermore, their results were normalized 
in terms of the gluon mean free path A. That's a problem for small coupling expansions because 
A is zero in (hard-thermal-loop resummed) perturbation theory! A popular model for A, often 
used as an example by BDMPS, is perturbative Coulomb scattering cut off by Debye screening, 
with the form 



Preprint submitted to Nuclear Physics A December 22, 2009 



where q ± is the transverse momentum transfer, n~T 3 is the density of plasma particles, and 
mu~gT is the Debye mass. The Debye effect screens electric fields, but it does not screen nearly- 
static components of magnetic fields. As a result, magnetic scattering can take place at much 
smaller q ± , and the actual formula representing perturbative scattering with screening effects is 



This is infrared log divergent, giving A = 0. 

The result I have to report J^t] is that, with minor tweaking of the way BDMPS-Zakharov 
write their formulas, one can eliminate all the above issues and then easily show that BDMPS- 
Zakharov reproduces AMY in the uniform, infinite medium limit. So, there is now a BDMPS- 
Zakharov formalism that handles non-static scatterers, and equivalently there is now a version of 
AMY for non-uniform media. See Ref. ]5|] (section II. B and the appendix) for details. The basic 
tweak to BDMPS-Zakharov is to write their formulas more generally in terms of the rate T e i for 
elastic scattering off of the medium instead of as density n times the corresponding cross-section 
cr, and to avoid normalizing quantities by A. 

2. Large logarithm approximation 

To fully evaluate the bremsstrahlung rate to leading order in a s requires somewhat compli- 
cated numerics, even in the case of uniform media. There is a great deal of simplification if 
one makes the additional approximation that the logarithm ln{E/T) of energy is large [while 
still treating a s ln(£'/r) <K 1]. With this large logarithm approximation, it is possible to get 
analytic results. To make practical use of a large logarithm expansion, however, one needs a 
next-to-leading log (NLL) result. Leading log cannot tell the difference, for example, between 



But these two logarithms are very different for realistic values of E/T. 

But, even if one has a NLL result, is it any use in practical situations, or is log(E/T) never 
large enough? A test is shown in Fig.Q] which gives results from Ref. |0]. The solid lines are 
results from a full leading-order evaluation of the g—>gg bremsstrahlung rate (equivalent to the 
sort of calculations first performed by Jeon and Moore 01 based on the AMY formalism), as 
a function of the bremsstrahlung gluon energy at = xE divided by T. The dashed lines show 
the result of the NLL calculation of Ref. [6]. The horizontal axes goes up to the very large, 
absurdly unrealistic value u/T = 10 5 just to verify that the curves do approach each other as 
ln(E/T) — > oo (for fixed x), as they should. The conclusion to take away from this plot is that, in 
the context of a weakly-coupled plasma, the large logarithm expansion to NLL is good to < 20% 
for lo > 1QT, which is much better than one might have feared. For u> ~ T, the NLL results is 
off from the full small-coupling result by roughly a factor of two. 

3. q in weakly-coupled plasmas 

Finally, I want to wrap up this potpourri by discussing the value of the jet broadening param- 
eter q in the limit of weak coupling, q is defined as the averaged squared transverse momentum 
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Figure 1: g— >gg bremsstrahlung rate in an infinite, uniform medium. Solid line is a full calculation to leading order in 
a s . Making the further approximation that ln(xE/T) is large, and working to NLL order, gives the dashed line. See Ref. 
((J for normalization of vertical axis. 




Figure 2: Elastic scattering of a high-energy particle (p) off of a plasma particle (pi). 



transfer Q\ per unit length to a high-energy particle traversing the plasma, so that Q\ = qL. 
It's of relevance to bremsstrahlung because the formation time depends on the collinearity of the 
bremsstrahlung gluon with the particle that emits it, and the degree of collinearity in turn depends 
on how much the particles are randomly deflected during the bremsstrahlung process. 
The squared transverse momentum transfer per unit length is simply 

q = J d 2 q ± -j^- q 2 ± , (4) 

where r e i is the rate for elastic scattering from plasma particles, as in Fig. [2] and q ± is the 
transverse momentum transfer in a single collision. When evaluated at leading-order in a s (and 
ignoring the running of a s ), this formula is UV log divergent. For the bremsstrahlung problem, 
however, what one actually needs is the result q{A) obtained by introducing a UV cut-off A on 
the q ± integration |H EE • 

The calculation of the dF e i/d 2 q ± needed in ©, based on the leading-order process of Fig. [2] 
has the form 

Jpi ~ j dq z j d 3 P2 §^/(/>2)[i ± /(pa - q)l (5) 

where f(p) is the Bose or Fermi distribution for finding a particle with momentum p in the 
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plasma. The 1 ± / factor is a Bose enhancement or Pauli blocking factor for the final state of 
the plasma particle. Most perturbative calculations of q have made additional simplifications to 
this formula. Some implicitly ignore the final-state 1 +/ factor, which is ignorable only when 
the dominant momentum transfers q are large compared to T, and even then only to leading-log 
order. Others replace 1 + f(j>2 - q) by 1 + f(pi), valid if the dominant q are small compared to 
T. To my knowledge, the result has only recently been evaluated using the full form © [9, [H. 
For A » T [but a s ln(A/T) <k 1], an analytic small-coupling result for UV-regulated q is given 
in Ref. [9]. As an example, for a purely gluonic plasma, it i^J 
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where 



c = 2exp(± -ye), 



ln[(*-l)!] 

k 3 



0.386043817389949. 



(6b) 
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and mu = gT is the Debye mass. The constant <x + can be related to certain generalizations of the 
Riemann £ function. I think the expression (0 for q is fun and interesting. But if you are of a 
more practical bent, you could have instead just done the integral (0 numerically. 

A presentation of q to leading order in <? s requires an important warning: Corrections which 
are formally higher-order in coupling, of order mu/T = 0(g), have been analyzed by Caron-Huot 
ifioll and are of order 100% for realistic couplings. 
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1 For comparison, if one made the q <K T approximation of replacing 1 + f(j>2 — 9) by 1 + f(pi) in (5), the result 
would be q g (A) = [f(2)ln^] 9s w i ■ If one instead made the q 2> T approximation of dropping the 1+/ term in j5J, 
the result would be §g(A) = [^(3) In ^1 ' >g J ■ A rough approximation one sees in the literature is to (i) replace the 
final factor of q\ in {4) by m^, which artificially eliminates the UV log divergence and gives q ^ tn^/A, and (ii) use 
the model form for the gluon mean free path A. This gives q g =* [ 3 f (3)j 9s ^ in the purely gluonic case, for which 
# = 9^(3)/2tt 4 in (TJ. 
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